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Abstract 

A series of short time stochastic resonance (SR) phenomena, realized in a bistable 
receiver, can be utilized to convey train of information represented by frequency-shift 
keying (FSK) signals. It is demonstrated that the SR regions of the input noise inten- 
sity are adjacent for input periodic signals that differ in frequency appropriately. This 
establishes the possibility of decomposing M-ary FSK signals in bistable receivers. Fur- 
thermore, the mechanism of the M-ary FSK signal detection via short time SR effects 
is explicated in terms of the receiver response speed. The short time SR phenomenon 
might be of interest for neuronal information processing in non-stationary noisy envi- 
ronments, regardless of the short timescale or the frequency jitter of stimulus. 

1 Introduction 

Stochastic resonance (SR) is now a well established phenomenon wherein the response of a 
nonlinear system to a subthreshold periodic input signal can be enhanced by the assistance 
of noise [1—7]. Since a single- frequency sinusoidal input conveys little information content, 
this effect has been extended to aperiodic (i.e. broadband) input signals, leading to the term: 
aperiodic stochastic resonance (ASR) [8-16]. From the point of view of information transmis- 
sion, the aperiodic information-bearing signal might be associated with analog (amplitude 
and frequency) modulated signals [17-21], or digitally modulated signals [10,13,21-25] within 
the context of SR effects. The digital pulse amplitude modulated signal is intensively inves- 
tigated for revealing new ASR phenomena or novel applications [10,13,22-25]. Recently, a 
new type of electronic receiver based on SR properties has been proposed for retrieval of a 
subthreshold frequency-shift keying (FSK) digitally modulated signals [19,21]. 

In the present paper, the input information sequences are also represented by equal-energy 
orthogonal signal waveforms that differ in frequency, i.e. the FSK signals [19,21,30]. But, the 
prototype SR model, i.e. an overdamped bistable system [1-6,9-11,13,18,22-24,31,32], is 
adopted as a nonlinear receiver that decomposes the received signals. Usually, conventional 
SR or ASR characterizes the SR phenomenon with a statistical measurement, resulting from a 
long-term observational data [1-16] . In contrast, we are more interested in the noise-enhanced 
effects occurring in each short-term duration of each symbol interval as the noise intensity 
increases, what we call the short time SR phenomenon in this paper. The short time SR effect 
is consistent with detecting weak signals from a short data record [31], storing information in a 
short-term memory device [32] and exploring transient stimulus-locked coordinated dynamics 
in terms of mechanisms of short-term adaptation in sensory processing [33] . From a series 
of short time SR effects, the input information contents can be deciphered at the output 
of the bistable receiver in terms of different signal frequencies. The SR effects realized 
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Figure 1: Bistable receiver and demodulator for M-ary FSK signals. 



in a nonlinear bistable receiver, as will be shown, are not sensitive to adjacent periodic 
signals with an appropriate frequency separation. In other words, the values of the noise 
intensity at resonance regions or points are close for two adjacent periodic signals. This 
nonlinear characteristic of a bistable receiver is then studied for transmitting M-ary FSK 
digitally modulated signals in detail. Furthermore, the response speed of bistable receivers, 
independent of the frequency of input signals, is theoretically deduced in Sec.|3 In view of 
the receiver response speed, the mechanism of the M-ary FSK signal detection is explained. 
Finally, we emphasize that the short time SR is not a trivial effect in signal detection, for 
example, in the detection of bipolar pulse signals with an unknown arrival time [27]. The 
short time SR effect may also be of interest in neurophysiology, even the stimulus exists in 
a short-term timescale or has a frequency jitter. Neuronal noise, in non-stationary noisy 
environments, might generically enhance the prediction of the arrival information via the 
short time SR effect, before the adaptive capabilities of neurons lowers their thresholds to 
maximize information transmission [28,29]. 

2 Bistable receiver and M-ary FSK signal transmission 

An input information-bearing sequence {J} is mapped onto the M-ary FSK signal S(t) as 

S(t) = Acoa(2irf m t),m=l,2, ...,M, (1) 

for (n — 1)T < f < nT, n = 1, 2, • • •. Here, A is the amplitude, {f m , m = 1,2, ... , M} denotes 
the set of M possible carrier frequencies corresponding to M — 2 k possible fe-bit symbols, 
and T is the symbol interval. The received signal R(t) is 

R(t) = Acos(2ir f m t + m ) + 77(f), (2) 

where <f) m are the phase shifts of carrier frequencies f m induced by the channel, and the 
background noise 77(f) is additive Gaussian white noise with autocorrelation (77^)77(0)} = 
2DS(t) and zero-mean. Here, D denotes the noise intensity. Next, R(t), as shown in Fig.^ 
is applied to a bistable dynamic receiver given as 

T a ^ = x(t)-^+R(t), (3) 

with receiver parameters r a > and X^ > [22]. Here, r a is related to the system relaxation 
time. The dynamics of Eq. J3J is derived from the symmetrical double- well potential Vq(x) = 
—x 2 /2 + x 4 /(4X£), having the two minima Vb(±X{,) = — Xjj/4. Parameters r a and X^ have 
the units of time and signal amplitude respectively, and define natural scales associated to 
the process of Eq. © [23] . 

We are interested in recovering the successive input information bits, from the observation 
of the system state x(t). Moreover, our focus is on the signal frequency rather than the wave- 
forms at the output of bistable receiver, due to the information contents being represented 
by different frequencies. In this paper, we numerically integrate the stochastic differential 
equation of Eq. @ using a Euler-Maruyama discretization method with a small sampling 
time step Af <C t q [37]. The demodulation method, as shown in Fig.^ is adopted with the 
zero crossing times N m . If the bistable receiver follows the subthreshold periodic signal cor- 
rectly by the assistance of noise, the zero crossing times N m of modulated signal 5(f) should 



2 




1 2 3 4 5 




1 2 3 4 5 

time (s) 



Figure 2: Time evolution of binary FSK signal transmission in a bistable receiver. (A) 
S(t) = 50cos(87rf) representing binary digit 1 (or 0). S(t), plotted by solid and dotted curves 
alternately in each symbol interval of T = 1 s, should be considered as segment signals rather 
than a continuous one in the whole transmission time. (B) The received signal R(t). Here, 
the phase shift (j> m induced by the channel is tt/6 and the noise intensity D = 0.8 V 2 /Hz. 
(C) The bistable receiver output signal x(t) with parameters r a = 1/3000 s and Xi, — 150 V. 
The decoded binary bits are depicted in terms of the zero crossing times. The sampling time 
At = 10" 3 s. 



be 2f m T in each symbol interval T. At the output of the bistable receiver, however, the zero 
crossing times N m will be in the vicinity of 2f m T even in the resonance region of noise, as 
shown in Fig. [21(C). For simplicity, the bytes represented by the corresponding input signals 
with frequency f m are decoded as 

tfm-l + fm)T < N m < (f m + f m+ l)T, 

or (2/ m - A/)T < N m < (2/ m + A/)T, (4) 

with the frequency separation A/ = f m +i — fm for m = 1, 2, . . .. Then, the output informa- 
tion sequence Y is decoded. Now, this system of Eq. J3J) with input digits and output digital 
readings, can be viewed as an information channel transmitting digital data. By comparing 
sequences / and Y, the measure of the percentage P of bytes correctly decoded will be used 
to quantify the performance of this nonlinear information channel. We shall show that this 
transmission of information can be assisted by additive noise addition in each symbol interval, 
a property we interpret as a short time SR effect. 

In Fig. tuning the noise intensity D so that the switching between wells is made to 
agree closely with the input periodic signal, resulting from the combined action of noise 
and periodic signal in a bistable receiver. This synchronization phenomena in each symbol 
interval of T is called the short time SR effect. In Fig. 0(A), the input periodic signals in 
each short time data record of T have the same frequency, which can represent a particular 
input binary sequence I = [1 1 • • • 1] or / = [0 • ■ • 0]. The resonance curves, illustrated in 
Fig. are distinct from the previous conventional SR form: the SR effects is explored in each 
symbol interval T, rather than the entire transmission time, recording the binary information. 
Numerical results of Fig. show that the percentage P of correct switching events, as the 
noise intensity D increases, presents a typical SR characteristic. Especially worthy of note is 
that the resonant regions or points of two adjacent frequencies are close. This indicates that 
the bistable receiver is not sensitive to the frequency jitter of input periodic signals. Thus, 
it is possible to use a bistable receiver to decode binary or 4-ary information represented 
by periodic signals with different frequencies, i.e. FSK digitally modulated signals. This 
possibility is immediately demonstrated in Figs. 0] and [5] with numerical experiments. In 
different symbol intervals of T, the bistable receiver can follow the input periodic signals 
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Figure 3: The percentage P of correctly decoded bytes as a function of noise intensity D. In 
each symbol interval of T, the input signal S(t) — 50 cos(27r f m t) has frequency f m (as marked 
in legends). The statistical values of P are computed numerically from 5000 transmitted 
codes. The bistable receiver is with parameters r a = 1/3000 s and X b = 150 V. The symbol 
interval T = 1 s and the sampling time At = 10~ 3 s. 




Figure 4: Plots of the percentage P as a function of D for the binary FSK signal transmission 
in bistable receivers. The statistical values of P are obtained from 5000 transmitted codes. 
The corresponding carrier frequencies are given in plots. (A) The receiver parameters r a = 
1/3000 s and X b = 150 V. T = 1 s and S(t) = 50 cos(2ir f m t). (B) The receiver are with 
parameters t q = 1/10000 s and X b = 380 V. T = 0.1 s and S(t) = 140 cos(2tt f m t) . The 
sampling time At = 10~ 5 s. 
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Figure 5: Plots of the percentage P as a function of D for the 4-ary FSK signal transmission 
in bistable receivers. The corresponding carrier frequencies are given in plots. (A) and (B) 
are with the same parameters as in Fig. respectively. 
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Figure 6: A sinusoidal signal with the zero-hold sampling. 




with adjacent frequencies in an appropriate noise intensity region. With the decision rule of 
Eq. PJ l. the percentage P is a non- monotone function of noise intensity D, resulting from a 
series of short time SR effects. For M > 4, the carrier frequency distance increases and a 
bistable receiver is not sufficient for decoding the received FSK signals. Naturally, a parallel 
bank of bistable receivers can be designed for this complicated task, with different receiver 
parameters and improved deciphered scheme. 

3 The mechanism of the FSK signal detection in a 
bistable receiver 

In this section we shall now attempt to explore the physical mechanism of the FSK signal 
detection in a bistable receiver, related to a an approximation of the nonstationary probabil- 
ity density of Eq. (j3J and its temporal relaxation. The temporal relaxation of nonstationary 
probability density, termed the receiver response speed Ai, will be demonstrated being in- 
dependence of the input signal frequency. The study of system response speed allows us to 
explicitly have a deeper understanding of the FSK signal transmission in a nonlinear bistable 
receiver. 

In numerical simulations, we numerically sample a sinusoidal signal with zero- hold and 
sampling time At <C r a [10,22-24,27]. Hence, in a sampling interval of At, the system of 
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Figure 7: The system response speed Ai (Hz) as a function of signal amplitude A (V) and 
noise intensity D (V 2 /Hz) for the receiver with parameters (A) r a = 1/3000 s and Xj, — 150 V 
and (B) parameters r a = 1/10000 s and X b = 380 V. 



Eq. is subjected to a constant amplitude s(t) = Ai (iAt < t < (i + l)Ai, i = 0, 1, 2, • ■ •), 
as shown in Fig. EI I n the presence of noise r](t), the statistically equivalent description for 
the corresponding probability density p(x, t) is governed by the Fokker-Planck equation 

dp(x,t) r d T , D d 2 n , 

where V'{x) = -(x - x 3 /Xj; + Ai) and the Fokker-Planck operator L F p = -§^V (x) + ^--^ . 
p(x, t) obeys the natural boundary conditions that it vanishes at large x for any t [40]. The 
steady-state solution of Eq. (JSJ , for a constant input at Ai , is given by 

p{x) = lim p{x,t) = Cexp[-T a V(x)/D], (6) 

t— >oo 

where C is the normalization constant [40]. As the sampling amplitude Ai varies, we en- 
counter the nonstationary solution p(x,t) of the Fokker-Planck equation, i.e. Eq. J^J. This 
analysis is performed in Appendix A. We show in Appendix A that the nonstationary so- 
lution p(x,t) can be expanded as an asymptotic representation by eigenfunctions Ui(x) and 
eigenvalues 

n 

p{x,t) =J2CiUi(x)exp[-T a V(x)/{2D)} exp[-Ait], (7) 

where Ci are normalization constants for i = 0, 1, 2, . . .. 

Specifically, Ai, termed the receiver response speed, is a measure of the slowest time 
taken by the bistable receiver to tend to the steady-state solution of Eq. ©. In other words, 
Ai is the speed of bistable receivers tracing the variety of input signals, whence our term 
"receiver response speed." Figure shows the behavior of the receiver response speed in 
the related regions. Note that the receiver response speed Ai, for a fixed bistable receiver 
with parameters r Q and A"b, is a monotonically increasing function of signal amplitude A and 
noise intensity D, but independent of the input signal frequency f m , as indicated in Eq. (|15f) . 
When the noise intensity D is too small, the receiver cannot follow the input signal correctly, 
as shown in Fig. [S] (B). However, in the resonance regions of noise intensity D, Ai is large 
enough to make the receiver output x(t) reach the steady-state, whereas the signal amplitude 
continually changes from ±A to =pA at different but adjacent carrier frequencies f m . This 
indicates detecting weak FSK modulated signals is possible in a bistable receiver, as seen in 
Fig- El (C). After the noise intensity D is beyond the resonance region, Ai is too fast to catch 
up with the fluctuations induced by noise, resulting in the loss of synchronization, as shown 
in Fig. El (D). 
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Figure 8: (A) A 4-ary FSK signal S(t) = 140 cos(2tt f m t) with carrier frequencies 40 Hz, 
60 Hz, 80 Hz and 100 Hz, representing 4-ary code words 00, 01, 10 and 11, respectively. The 
receiver with parameters r a = 10~ 4 s and Xi, = 380 V. (B), (C) and (D) are plots of receiver 
outputs x(t) at D = 0.25 V 2 /Hz, D = 1.0 V 2 /Hz and D = 2.4 V 2 /Hz. The erroneous 
decoded bytes are bracketed. Here, T — 0.1 s and At = 10~ 5 s. 



The receiver response speed Ai contributes to the mechanism of transmitting weak FSK 
modulated signals in a bistable receiver: Within a reasonable region of noise intensity, the 
receiver responds synchronically to the input periodic signal, regardless of the short-term 
timescale of symbol interval T or the frequency jitter of signals (i.e. the different carrier 
frequencies f m ). Therefore, the SR effect realized in a symbol duration T, i.e. the short time 
SR, can be utilized to convey or store the digital data [31,32]. Additionally, noise-enhanced 
frequency discrimination was reported in Ref. [26], and the frequency robust characteristic 
of bistable model is also and verified in Ref. [27] in detail. 



4 Discussion 

A detection strategy of subthreshold FSK digitally modulated signals in a nonlinear bistable 
receiver is analyzed. It was numerically demonstrated that the resonance values of noise 
intensity appear closely for FSK modulated signals with adjacent frequencies. Assuming a 
demodulation method of the zero crossing times, a series of SR effects appearing in each 
symbol interval, what we call the short time SR, provide the possibility of detecting sub- 
threshold M-ary FSK signals in bistable receivers. In order to understand the mechanism 
of the FSK signal detection more deeply, we introduce the receiver response speed, i.e. a 
theoretical measure for the receiver tending to a steady-state. In the resonance region, the 
receiver response speed is fast enough to trace the variety of input periodic signals, regardless 
of the frequency difference and the short-term timescale of symbol interval. 

Finally, we argue that the short time SR effect might be an optional detection strategy in 
neurodynamics. Neural noise will play a positive beneficial role for information processing, in 
the case of the stimulus exists in a short time duration or has a jitter in frequency [26,31,32]. 
Even under the assumption that neurons have sufficient adaptability in sensory systems 
[28,29], neural noise might provide an enhancement via the short-time SR phenomenon 
before the neuron adaptability acts. This is an open question and currently under study. 
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A Receiver response speed and nonstationary probabil- 
ity density model 

In Eq. (J5J, the Fokker-Planck operator Lpp = gjV'(x) + ^--§^s is not a Hermitian operator 
[40]. We rescale the variables as 



X b = X b /y/D/r a , A, = Ai/y/D/ra, 

t = t/ra, y = x/^D/ Ta , (8) 

Eq. JS} becomes 

where V'(y) = — (y — y 3 /X 2 + Ai). The steady-state solution of Eq. © is given by 

p{y) = lim p(y, T ) = Cexp[-V(y)], (10) 

where C is the normalization constant. A separation ansatz for p(y,r) [40], 

p(y, t) = u(y) exp[-V(y)/2] exp(-Ar), (11) 

leads to 

Lu = -Ait, (12) 

with a Hermitian operator L = Jp- — [jV' 2 (y) — \V"(y)\. The functions u{y) are eigenfunc- 
tions of the operator L with the eigenvalues A. Multiplying both sides of Eq. Ijl2|l by u(y) 
and integrating it, yields 

, _ J^W 2 (y) + uHyW 2 (y) - \V''{y)]}d y 

J-oo U \V) d V 

where eigenfunctions u(y) satisfy the boundary conditions of u(y) and u'{y) tending to zero 
as y — > ±oo. The eigenvalue problem of Eq. (|12[) is then equivalent to the variational 
problem consisting in finding the extremal values of the right side of Eq. (|13|l [24,40]. The 
minimum of this expression is then the lowest eigenvalue Ao = 0, corresponding to the steady- 
state solution of Eq. (fTT)fl [40]. We adopt here eigenfunctions u(y) = p(y) exp[— V(y)/2] and 
p(y) ^ 0, Eq. {E3> becomes 

r +0 ° i 

X = { J co { p,2 ^ + 2 p2{y)V ' 2{y) - 

\[V'(y)p 2 (y)]'}eM-V{y)]dy} 

/+OC 
P 2 (y)eM~V(y)]dy}. (14) 
CO 



Since 



[V'(y)p 2 (y)}'exp[-V(y)]dy 

O 

= V'(y)p 2 (y)eM~V(y)}\±™ 

/ + CO 
p 2 (y)V' 2 (y)eiq ? {-V(y)}dy 
-CO 

f + oo 

p 2 (y)V' 2 (y)exp[-V(y)}dy, 
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Eq. lfH| can be rewritten as 



A = — — . (15) 

Loo P 2 (y) exp[-V(y)]dy 

Assume p(y) = do + d\y + . . . + d n y n and the order n is an integer, we obtain 

([K]-X[M]){d} = 0, (16) 

with eigenvectors {d 1 } = [d l ,d\, . . . d l n ] corresponding to eigenvalues {A} = [Ao, Ai, . . . , A„] 
for i = 0, 1, .. .n. The integer n is not increased in the iterative process until the preceding 
values of A^ approximate the next ones within the tolerance error. The elements of matrices 
[M] and [K] are 



/+oo 
y i+ Uxp[-V(y)}dy>0, 
-OO 



-\-oc 



kij = / ijy i+j - 2 eM-V(y)]dy > 0, 



where i, j = 0, 1, . . . n. The matrix [M] is positive definite and the matrix [K] is semi-positive 
definite. The minimal eigenvalue Ao is zero. The minimal positive eigenvalue Ai describes 
the main speed of the system tending to the steady state solution of Eq. i|10fl , what we call 
the receiver response speed. 

From Eq. (|16|l . we can obtain the eigenfunctions Ui(y) = pt(y) exp[— V(y)/2] correspond- 
ing to the eigenvalue Ai for i = 0, 1, . . . n, where Pi{y) = d l +d\y+. . .+d l n y n . The eigenvectors 
{d 1 } = [d l Q ,d\, . . . d l n ] are normalized. Because L is a Hermitian operator, eigenfunctions Ui{y) 
and Uj (y) are orthogonal 

Ui{y)uj{y)dy = Sij, (17) 



where i,j = 0, 1, . . . ,n. Hence, p(y,T) can be expanded, according to eigenfunctions Ui(y) 
and eigenvalues A^, as 



Piy, T ) =^2 CiUi(y) exp[-V(y)/2] exp[-Ajr], 



i=0 



where Ci are normalization constants deduced from the orthogonal condition of eigenfunctions 
[40]. Note the scale transformation in Eq. JHJ, p{x, t) can be represented as 

n 

p(x, t)=J2 CiUi(x) exp[-T a V(x)/(2D)] exp[-A,t], (18) 

i=0 

with the real eigenvalues A; = A^/r a in the timescale of t and derived from Eq. ffify . 
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